Effect of non-linear interface kinetics on coarsening phenomena 
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Coarsening kinetics is usually described using a linear gradient approximation for the underlying 
interface migration (IM) rates, wherein the migration fluxes at the interfaces vary linearly with the 
driving force. Recent experimental studies have shown that coarsening of nanocrystalline interface 
microstructures is unexpectedly stable compared to conventional parabolic coarsening kinetics. Here, 
we show that during early stage coarsening of these microstructures, IM rates can develop a non- 
linear dependence on the driving force, the mean interface curvature. We derive the modified mean 
field law for coarsening kinetics. Molecular dynamics simulations of individual grain boundaries 
reveal a sub-linear curvature dependence of IM rates, suggesting an intrinsic origin for the slow 
coarsening kinetics observed in poly crystalline metals. 

PACS numbers: 81.07.Bc, 81.10.Aj, 61.72.Cc, 61.82.Rx, 61.72.Mm, 68.35.-p 
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Material properties of most inorganic polycrystals de- 
pend on the underlying interfacial microstructure, in par- 
ticular the grain size. The positive free energy of the in- 
terfaces provides a universal driving force for grain coars- 
ening such that a curved interface moves towards its cen- 
ter of curvature in order to decrease the total system 
energy. A quantitative understanding of the curvature 
dependence of the interface migration (IM) rate is cen- 
tral to much of material processing as it determines the 
overall microstructure coarsening kinetics. 

In systems where interface motion is activated, the 
atomic-exchange across the interface determines the IM 
rate, v. Absolute reaction rate theory based on such 
atomic hopping events yields the dependence of the IM 
rate on the driving force p [l[ , 
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where b is the interface displacement per event, N is the 
number of event sites, v a is a jump frequency character- 
istic of the underlying lattice and AG a = AQa — TASa 
is the activation barrier for each hopping event. System 
free energy change per event Ag is the energy provided 
by the driving pressure p to effect the atom-exchange, i.e. 
Ag = pu) m , where u> m is the activation volume associated 
with each event. At high temperatures, (3pu) m <C 1 and 
Eq. [T]can be linearized to express the IM rate in terms of 
the driving force, v f» Mp. The constant of proportion- 
ality M is the interface mobility, which is predominantly 
Arrhenius with temperature 
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and M = fibNv u)r> 
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During coarsening, the capillary driving force on each 
interface segment is the product of the interface stiff- 
ness r and its mean curvature k, or the weighted mean 
curvature k 7 = Tk Q. The IM rate now increases 
with the mean interface curvature [!, H]> HI- The cur- 
vature dependence implies that material systems with 
ultrafmc/nanocrystalline (nc) grain sizes possess an in- 
herently high driving force for coarsening. This is often 
an unwanted outcome during thermal annealing, as the 
qualitatively superior thermomechanical and transport 
properties of these microstructures are offset by their in- 
stability with respect to coarsening. 

However, a poorly understood feature of nanocrys- 
talline microstructures is that the coarsening is anoma- 
lously slower than expected, almost linear in time, before 
it transitions to the conventional parabolic growth 0, 0] . 
This behavior has been attributed to grain size depen- 
dent extrinsic effects on interface motion, such as en- 
hanced vacancy and triple junction drag @, S, E3], 
solute segregation and particle incidence. In this ar- 
ticle, we present an alternative framework for coarsening 
behavior based on intrinsically non-linear IM rates, and 
discuss the implications of this behavior for the stability 
of nanocrystalline microstructures. 

The rationale is the observation that capillary driving 
pressures during early stages of coarsening in nanocrys- 
talline microstructures are large enough such that the 
linearized rate theory begins to break down. This is evi- 
dent from Table HI a list of the critical driving forces p cr 
in various polycrystalline metals at which [3p cr uj m = 0.1. 
The definition corresponds to a processing temperature 
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T = 0.7 T m and assumes single-atom hops across the 
interface, uj rn = Q (T m is the bulk melting point and 
O is the atomic volume). Since the atomic-scale mech- 
anism can in general involve more than one atom, as 
in correlated or military atom transfers across the inter- 
face 0, 12 1 , io m > n and the reported values of p cr are 
upper bounds. 



be generalized as 



A 



Metal system 


Volume ui m (A) a 


p cr (MPa) 


R cr (nm) 


Aluminum 


16.6 


54.4 


25.1 


Copper 


11.8 


87.7 


51.5 


Nickel 


10.9 


153.2 


18.3 


Lead 


30.3 


19.1 


146.4 


Gold 


17.1 


76.1 


36.8 


Silver 


16.9 


70.8 


39.6 



TABLE I: Critical driving forces for interface migration and 
critical grain radii in polycrystalline metals at which lin- 
earized reaction rate theory breaks down. The activation vol- 
ume for IM is assumed to be the atomic volume, iu m = fi. 



The critical driving force can also be used to dchne 
a critical grain size, R cr , as the average grain size in 
polycrystals is a measure of the capillary driving force 
per interface segment. The mean field relation takes the 
form re 7 = TB/R, where the topological parameter B 
captures the effect of the interface network [13J]. Table 1 
also lists the lower bound for the associated critical grain 
sizes below which the non-linearities become important. 
The interface stiffness is assumed to T = 0.5 J/m 2 , while 
the topological parameter is based on Potts model grain 
growth simulation-based comparison between shrinking 
of an embedded sphere (B = 47r) and a 3D polycrystal, 
B ~ 2.8 [ADR, to be published] . These critical grain sizes 
lie well within the range of those found in nanocrystalline 
microstructures in these metals, emphasizing the role of 
non-linear interface migration kinetics in determining the 
overall coarsening kinetics and therefore their thermo- 
mechanical stability. 

The extension of the reaction rate framework to non- 
linear IM is not straightforward. The activation barrier 
is now comparable to the driving force, the bias in the 
energy landscape (inset, Fig[T]). As a result, the nature 
of the bias becomes important. While several minimum 
energy paths (MEP) are possible depending on the effect 
of the driving force on the initial, activated and final 
states, we limit ourselves to two extreme scenarios : a) 
the forward activation barrier is unaltered while the final 
state is lowered by the applied bias, or b) the bias is 
distributed equally between the initial and final states 
such that the forward barrier is altered. In general, the 
energy of the activated state in the resultant MEP is also 
changed due to the driving force. Making the simplifying 
assumption that this energy change is small, Eq. [1] can 
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where A = /3M/u> m (181 ]. Equation [3] is similar in form to 
the Butler- Volmer equation in electrokinetics for transfer 
of charged species [14[ . 

The parameter £ captures the combined effect of the 
degree of symmetry in the energy landscape imposed by 
the driving force, and the change in energy of the acti- 
vated state. The two extreme cases correspond to £ = 
and £ = 0.5. The £ = scenario can arise if the driving 
force is distributed asymmetrically such that its effect is 
restricted to the final state. The scenarios < £ < 0.5 
result in a decrease in effective forward barrier with in- 
creasing driving force. For £ = 0.5, the bias is distributed 
symmetrically about initial and final states such that the 
forward activation barrier is AG a — AG a — (Ag)/2 (see 
inset Fig. [J) . Changes in the activated state can be read- 
ily absorbed into an equation of this form. For example, 
if the MEP for the £ = 0.5 scenario is modified such that 
the activated state also increases by an equal amount 
(Ag/2), we recover the £ = scenario. 




r=R/a 



FIG. 1: (color online). Normalized migration rate as a func- 
tion of reduced grain size f = R/a, for linear (circles), sub- 
linear (Eq.|4]with £ = 0, squares) and super- linear (Eq.|4]with 
£ = 0.5, diamonds) driving force dependence of IM rates, (in- 
set) Schematic of the energy landscapes for in the absence 
and presence of a driving force, for the two values of f . 

Limiting our analysis to the two scenarios, the IM rates 
can be expressed as 

v(0) = A(l - e -P pUm ) and v(0.5) = 2A sinh(-/3pw m /2). 

In the limit f3puj m 3> 1, v becomes independent of £, 
justifying the linearization (Eq.[T]) at small driving forces. 
In the mean field limit, the IM rate depends on the grain 
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size, v = dR/dt and we arrive at the relation between the 
form of the energy landscape and the grain size evolution, 



dr 

a— = A 
dt 



(4) 



Here f = R/a is the dimensionless grain size. The 
normalization factor a = BY(3uj m is a fundamental mi- 
crostructural length scale that is related to the critical 
grain size, R cr = 0.1a. Figure [1] shows the IM rates pre- 
dicted by Eq.[4]as a function of reduced grain size for the 
two scenarios. The linear approximation is also plotted 
for comparison. At small grain sizes such that f — > 1, the 
variation of IM rate is significantly non-linear with the 
driving force. The linearized relation is an overestimate 
(sub-linear IM rate) or an underestimate (super-linear 
IM rate) depending on the value of £, underscoring the 
role of the nature of the bias in the energy landscape. 

The non-linearities in IM rates will modify the over- 
all coarsening kinetics. For an isotropic interfacial mi- 
crostructure with an initial grain size fj, the mean field 
governing equation for the final grain size is 
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Putting y = 1/r, with Y f{ ^ = (l/r /(?) ) and Y. L = (1/r,), 

(5) 
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This form is convenient for the following standard 
Maclaurin expansion (Ref. p. 804): 
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Here, B m (X) are the Bernoulli polynomials. They are 
tabulated (Ref. [15], p. 809); for example, 



B {\) = 1, Si (A) = — i + A and B 2 (A) 



A + A 



Also, B m {\ - A) = {-\) m B m (\). Thus, the integral in 
Eq. [5] can be written as 

m=0 JY f(Z) 

and the grain growth equation becomes 
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with C m (0 = 



(-l) m B m (Q 
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Remembering r—(l/y), we arrive at the coarsening law 
2At 



h(f) 
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At high temperatures and large grain sizes F 3> 1, Eq.Q] 
can be linearized and we recover a linear relation between 
IM rates and the driving force, dr/dt w A/af. The re- 
sultant governing equation for coarsening kinetics is, 
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rdr = r f(p) - r\ 
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where is the final grain size. As expected, when 

ff( v ) fi, Ff(p) Vt- Combining Eqs. [6] and [7] yields 
the difference in grain size due to non-linear IM rates, 
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Coarsening kinetics described by Eq. [6] is shown in 
Fig. [21 a plot of fj(o) and Fjm.s) against 2A(t/a) — 
h[ff(£)] — h[l], for ^ = 1. Contribution of 0(f -2 ) and 
higher terms is negligible and ignored. Grain size pre- 
dicted by parabolic coarsening kinetics due to linear IM 
rates, r/( p ), is also shown for comparison. For sub-linear 
IM rates (£ = 0), the linear term dominates during early 
stage coarsening. The negative deviation from classi- 
cal parabolic coarsening kinetics increases linearly with 
grain size. We approach parabolic coarsening kinetics 
for f/(o) S> 1) yet the final grain size is smaller com- 
pared to classical parabolic coarsening - the coarsening 
is suppressed. The effect is exaggerated at smaller initial 
grain sizes (inset B in Fig. [5]); coarsening kinetics is in- 
creasingly linear (f/(o) < !)■ For super-linear IM rates 
(£ = 0.5), the positive deviation from parabolic coars- 
ening increases logarithmically and is much slower; the 
increase is substantial only for decades increase in grain 
sizes. The deviation is enhanced for fi < 1, as IM rates 
are super-linear for larger range of grain sizes. This can 
be seen in inset B in Fig. [21 for fi = 0.1. 

Our analysis shows that a fundamental understand- 
ing of overdriven interface motion is critical for predict- 
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FIG. 2: (color online). Reduced grain size f/ vs. normalized 
time 2A(t/a) predicted by Eq. [6] for initial grain size f; = 
1. Circles, squares and diamonds represent linear, sub-linear 
(£ = 0) and super-linear (£ = 0.5) IM rates-based coarsening, 
respectively. Inset A shows an enlarged view for 1 < f; < 2. 
Inset B shows the variation for fi —0.1 and 0.1 < rj < 1. 



ing coarsening kinetics in nanocrystalline interfacial mi- 
crostructures. Recent molecular dynamics (MD) simula- 
tions of flat bicrystals in pure Al have been performed at 
small and large driving forces, offering a basis for under- 
standing the coarsening kinetics of grain boundary mi- 
crostructures. The mobility of a flat = 38.2° <111> tilt 
misorientation grain boundary has been extracted in the 
zero driving force limit, and also under the influence of a 
bulk body force. Both studies were performed using an 
embedded-atom-method (EAM) framework for the inter- 
atomic potentials, justifying the comparison. The former 
is based on the random walk of the mean grain boundary 
position due to the uncorrelated thermal noise in the sys- 
tem [r| and was extracted at a temperature T — 750°K, 
the latter on a synthetic driving force due to an orien- 
tation dependent bulk energy term at T 



800°K [17 1 • 

The mobility of the synthetically driven grain bound- 
ary was extracted using a driving force 0.025 eV/atom. 
Using a conservative estimate of the activation volume, 
Lu, m = QaIi the driving force is still well past our critical 
value, 0pcj m — 0.4. Furthermore, the driving force for 
this boundary is of the order of the activation energy for 
migration of this boundary (~ 0.02 eV), extracted in the 
zero driving force limit [ZTT and MU, to be published]. 

Comparison of the absolute mobilities grain boundaries 
extracted using these two techniques reveals that zero 
driving force limit mobility is faster by almost an order 
of magnitude (4.4 x 10~ 7 vs 6.0 x 10 -8 m 4 J _1 s -1 )]), even 
though it was extracted at a slightly lower temperature. 
While systematic studies are necessary to understand the 
rather large difference (changes in energy of the activated 



state or large activation volume for migration), the de- 
crease in mobility and therefore the IM rates due to non- 
linearities suggests that coarsening of nanocrystalline is 
intrinsically suppressed, also confirmed in experiments. 
This intrinsic effect must be factored in before ascrib- 
ing the slow coarsening rates of these microstructures 
to extrinsic effects on interface motion. Meso-scale grain 
growth simulations are necessary to understand the com- 
bined effect. 
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Change in the energy of the activated state can also be 
absorbed in a similar form. If this increase is Cpo> mi the 
pre-factor A for the IM rate is scaled by e ~ (l3p " m . 



